A spectral Galerkin method for the the coupled Orr–Sommerfeld and induction equations for free-surface MHD

We develop and test spectral Galerkin schemes to solve the coupled Orr–Sommerfeld (OS) and induction equationsfor parallel, incompressible MHD in free-surface and fixed-boundary geometries. The schemes’ discrete bases consistof Legendre internal shape functions, supplemented with nodal shape functions for the weak imposition of thestress and insulating boundary conditions. The orthogonality properties of the basis polynomials solve the matrix-coefficient growth problem, and eigenvalue–eigenfunction pairs can be computed stably at spectral orders at leastas large as p = 3,000 with p-independent roundoff error. Accuracy is limited instead by roundoff sensitivity dueto non-normality of the stability operators at large hydrodynamic and/or magnetic Reynolds numbers (Re, Rm &4× 10). In problems with Hartmann velocity and magnetic-field profiles we employ suitable Gauss quadrature rulesto evaluate the associated exponentially weighted sesquilinear forms without error. An alternative approach, whichinvolves approximating the forms by means of Legendre–Gauss–Lobatto (LGL) quadrature at the 2p − 1 precisionlevel, is found to yield equal eigenvalues within roundoff error. As a consistency check, we compare modal growthrates to energy growth rates in nonlinear simulations and record relative discrepancy smaller than 10 for the leaststable mode in free-surface flow at Re = 3 × 10. Moreover, we confirm that the computed normal modes satisfyan energy conservation law for free-surface MHD with error smaller than 10. The critical Reynolds number infree-surface MHD is found to be sensitive to the magnetic Prandtl number Pm, even at the Pm = O(10) regimeof liquid metals.